Ionization branching ratio control with a resonance attosecond clock 
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We investigate the possibility to monitor the dynamics of autoionizing states in real-time and 
control the yields of different ionization channels in helium by simulating XUV-pump IR-probe 
experiments focused on the N=2 threshold. The XUV pulse creates a coherent superposition of 
doubly excited states which is found to decay by ejecting electrons in bursts. Prominent interference 
fringes in the photoelectron angular distribution of the 2s and 2p ionization channels are observed, 
along with significant out-of-phase quantum beats in the yields of the corresponding parent ions. 
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The evolution of valence electron wave packets in 
atoms, molecules and solids takes place on a timescale 
ranging from tens of attoseconds to few femtoseconds 
For example, the sudden removal of an electron in CO2 
and N2 molecules [3] initiates multielectron dynamics 
that unfolds on the attosecond timescale, a localized va- 
cancy propagates across the full length of a molecule as 
large as a tetrapeptide [4] within just ~ 1 fs, and a photo- 
electron escapes through the surface of solid tungsten [E| 
in 150 as or less. Recent advances in the generation of 
ultra-short pulses provide the tools necessary for a 
time-resolved pump-probe investigation of such dynam- 
ics and bear the promise of its control. 

Most of the reactive processes promoted by high elec- 
tronic excitation, like resonant multiphoton atomic ion- 
ization 0|, ultrafast electron-transfer [Toj , and molecu- 
lar dissociative photoionization [lll-[l3|, involve the for- 
mation of metastable, multiply excited states as a cru- 
cial intermediate step. These metastable states differ 
from bound excited states in that they can decay on a 
timescale that is comparable to the characteristic time 
evolution of the electronic wave packet itself. Their dy- 
namics is an essential ingredient of the rearrangement 
of correlated multielectron wavefunctions [14l-ll6j, and 
is thus of particular relevance for their eventual con- 
trol. In the present letter, we simulate a realistic XUV- 
pump IR-probe experiment focused on the N=2 ioniza- 
tion threshold of helium, the prototype of a multielectron 
system and the simplest neutral atom in which autoion- 
izing states arise. We show that the coherent superposi- 
tion of doubly excited states (DES) created by the XUV 
pulse ejects electrons in bursts rather than continuously, 
and we demonstrate that it is possible to exploit this dy- 
namics to effectively control the branching ratios of the 
different ionization channels. 

In our simulations, the time-dependent external field 
comprises an XUV-pump pulse followed by an intense IR- 
probe pulse, both with a Gaussian envelope. The XUV- 
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FIG. 1: Charge density after the XUV-pump pulse, at small 
(top row) and large (bottom row) radii. At each breathing cy- 
cle, the metastable wave packet, formed by a coherent super- 
position of doubly excited states, ejects a burst of electrons. 
The peak of the free electron density originating close to the 
nucleus results in a wavefront which propagates outward at 
almost constant speed, up to very large distances. 



pump pulse is 385 as long (full width at half maximum of 
the intensity), with the energy peaked at 60.69 eV, and an 
intensity of 2 • 10 13 W/cm 2 . The probe is a Ti: Sapphire 
800 nm (1.55 eV) pulse, 3.77 fs long (fwhm), with an 
intensity of 10 12 W/cm 2 . The XUV pulse populates a 
coherent superposition, \rfp), of 1 P° DES below the N=2 
threshold, mainly those belonging to the principal sp+ 
series [l7| ]: 

W P )-Y,\ S Pn)Cne-^ t / h (1) 

n 

where Z n = E n — iT n /2 is the complex energy of the 
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sp+ resonance, with position E n and width Y n . The lo- 
calized part of each term in this series is approximately 
represented by a symmetric linear combination of sp con- 
figurations, sp+ oc 2snp + 2pns [I?}- As a consequence, 
the localized part of \$p) is characterized by a symmet- 
ric breathing of p and s orbitals coupled to the 2s and 
2p parent ions, respectively: With the present pulse pa- 
rameters, the two lowest DES in the sp + series, sp~2 and 
sp^, which lie -5.04 eV and -1.69 eV below the N = 2 
threshold, with lifetimes of — 17.6 fs and — 80.3 fs, re- 
spectively, are by far the most populated ones. For sev- 
eral tens of femtoseconds, these two states dominate the 
dynamics of the met ast able wave packet. 

For the present simulation, the time-dependent 
Schrodinger equation is integrated numerically with an 
exponential propagator 

ij)(t + dt) = exp [-iH(t + dt/2)dt/h] ^(t), (2) 

where H(t) is the atomic Hamiltonian in velocity gauge. 
The wavefunction if; is expanded in a multi-channel close- 
coupling B-spline basis with total angular momentum up 
to L = 6 and the right hand side of Eq. (|2]) is evaluated 
with the Arnoldi algorithm. Each subspace with definite 
angular momentum L comprises the 1s</>l, 2s</>l, 2p0L+i, 
and 2pcj)L-i (for L > 0) close-coupling channels, where 
the notation nl(j)i> indicates that one electron is frozen 
in the nl He + orbital, while the other electron has the 
orbital angular momentum V . In the S symmetry, the 
basis also includes the Hartree-Fock ls 2 HF configuration 
for a better representation of the ground state. The ra- 
dial part of the atomic orbitals is expanded in a B-spline 
basis of order 10, with an asymptotic spacing between 
consecutive nodes of 0.5 Bohr radii, up to a given max- 
imum radius R. To compute the yield of the excited 
ions, a box with R — 400 Bohr radii was found to be 
sufficient, while for the partial differential photoelectron 
angular distributions (PDPAD) a larger box, R — 800 
Bohr radii, was used. In order to prevent reflections at 
the box boundaries, a channel-specific absorbing poten- 
tial V is included in the Hamiltonian: 

V = cJ2v a , V a =P a (r - R ) 2 9(r - Ro)P a , (3) 

a 

where the sum runs over all channels, P a is the projector 
onto the close-coupling channel a, 9{x) is the Heaviside 
step function, and c is a complex coefficient chosen as 
c = — (1 + 5z)10 -4 . The radius beyond which the po- 
tential is active, Ro, is set to — 100 Bohr radii from the 
box boundary. The absorbing potential V allows one 
to record the annihilation rate in each channel and to 
reconstruct the yields of all the parent ions. The photo- 
electron distribution in a channel, identified by a parent 
ion a = Is, 2s, 2p, is obtained by projecting the prop- 
agating wave function ^(t) onto the helium scattering 
states which satisfy incoming boundary conditions in all 
open channels but a [18|: 
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P a (E,Q)= 



(4) 



In Eq. (j4]), E and Cl denote the photoelectron energy and 
propagation direction, and the sum runs over the projec- 
tion m of the angular momentum of the electron in the 
parent ion, its spin cr, and the spin of the photoelectron 
c f . The scattering states are computed with the B-spline 
K-matrix method, a well-established configuration inter- 
action technique for the single ionization continuum [l9| . 

In FiglUi we show the electron density up to 700 Bohr 
radii at t=14.51 fs after the pump pulse. It consists 
of distinct wavefronts, spreading out with virtually con- 
stant speed, separated by time intervals which corre- 
spond closely to the beating period between the spj and 
the sp% resonances. In other words, the metastable wave 
packet decays by ejecting electrons in isolated bursts. 
This peculiar "cresting" behavior [20| can be understood 
in terms of interference between the long range part of 
the wave functions describing the decaying sp^ and sp% 
states. A more mechanistic interpretation, however, is 
possible. The Auger decay of DES is known to be trig- 
gered by electronic correlation; one of the electrons trans- 
fers part of its excitation energy to the other, which 
in turn is ejected into the continuum. Pisharody and 
Jones provided a spectacular and extreme example of 
this mechanism [2l[ ; they showed that the decay of some 
autoionizing states of helium, where both electrons are 
highly excited, takes place through a single violent e-e 
collision. A similar picture applies also when only one of 
the two electrons is highly excited [22]. In this case, the 
autoionization is found to take place at the encounter 
of the external electron satellite with the excited core. 
In the present case, though, neither of the two electrons 
is highly excited. In fact, the metastable wave packet 
has the smallest excitation possible, it lacks a clear semi- 
classical analogue, and the two electrons are constantly in 
close interaction. To investigate whether the collisional 
point of view still retain any validity, we traced the po- 
sition of 15 consecutive wavefronts in the time interval 
from 10 fs to 30 fs after the pump pulse, and extrapo- 
lated their evolution backwards in time to the moments 
at which they were created in the vicinity of the nucleus. 
The panels in the first and last columns in Fig. [T] cor- 
respond to two selected consecutive times at which a 
wavefront originates close to the nucleus, 14.51 fs and 
15.63 fs, while the central column corresponds to a time 
halfway between these two. In the upper row of Fig. [1] we 
show the electron density within 15 Bohr radii from the 
nucleus, which demonstrates its breathing motion. At 
t = 14.51 fs (a) the central part of the wave packet is at 
the peak of its contraction. At t = 15.09 fs (b) it reaches 
its maximal expansion. Finally, at t = 15.63 fs (c), it is 
contracted again. The instants at which the wavefronts 
are born in the vicinity of the nucleus therefore corre- 
spond closely to the stages of maximum contraction of 
the localized part of the metastable wave packet. This 
evidence supports the idea that the collisional descrip- 
tion of the autoionization dynamics of the DES of helium 
is indeed applicable down to the least excited ones. In 
the present case, though, it is not the encounter between 
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FIG. 2: Partial differential photo-electron spectra in the 2s 
(top row) and 2p (bottom row) ionization channels after the 
XUV-pump pulse (a,d) and after the IR pulse with two dif- 
ferent time delays At between pump and probe pulses sepa- 
rated by half the IR period: 15.53fs (b,e), and 16.87fs (c,f). 
x-axis: cosine of the photoelectron ejection angle with re- 
spect to the laser polarization; y-axis: photoelectron energy 
in atomic units. The interplay between the direct ionization 
by the XUV and the multiphoton ionization of the DES due 
to the IR-probe results in prominent interference fringes, with 
a characteristic energy spacing Ae = 27rh/ At. 

otherwise well-separated electrons [U, [22| that triggers 
the decay, but rather the squeezing of two electrons in 
constant close interaction. 

The XUV pulse has another major effect: it causes 
the sudden ejection of electrons in the Is as well as in 
the 2s and 2p channels. With the present choice of laser 
parameters, the proportion between the direct ionization 
in the Is channel, the population of DES and the direct 
ionization in the N=2 channels is roughly 10:1:0.01. In 
Fig. [2^ and Fig. [2]i we show the photoelectron angular 
distributions in the 2s and 2p channels, respectively, im- 
mediately following the XUV pulse, as functions of both 
the electron energy (y-axis) and the cosine of the angle 
between the electron propagation direction and the po- 
larization of the laser (x-axis). In the 2s channel, one 
recognizes the characteristic p distribution, proportional 
to cos 2 (the amplitude is odd), while in the 2p channel 
the angular distribution results from a combination of s 
and d waves (the amplitude is even). 

At the intensity considered, the IR-probe pulse has 
little effect on the ground state. It has a profound ef- 
fect, however, on the DES. The population of the 1 P° 
DES is partly redistributed among other DES with sev- 
eral different symmetries, and partly promoted to the 
continuum, mainly to the N = 2 channels. With an 
intensity of 10 12 W/cm 2 , the interaction of the system 
with the IR-probe pulse is a typical multiphoton pro- 
cess, where up to four IR photons are absorbed. As a 
consequence, the yield of the 2s and the 2p parent ions 
increases roughly by a factor of two, corresponding to 
~ 1% of the population of the DES. With more intense 
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FIG. 3: Yields of the 2s and 2p He + excited parent ions 
as functions of the time delay between the pump and probe 
pulses. Both yields are modulated by large quantum beats, 
due to the interplay between spj and sp% 1 P° DES, which are 
out of phase by as much as 60°. The continuous curves are 
obtained by fitting the computed points with a sine function 
plus a quadratic background. 

probe laser pulses, the yield of the excited He + ions can 
be substantially increased. IR laser pulses with a peak 
intensity of 10 13 W/cm 2 are routinely produced and pre- 
liminary calculations indicate that, at this intensity, the 
yield of N = 2 parent ions increases by more than one 
order of magnitude. The indirect multichannel ionization 
of DES is interesting because, by tracking the ionization 
yields in separate channels, one can follow the sharing of 
both the final energy and angular momentum between 
the two electrons, and hence the real-time evolution of 
electron-electron correlation in coherently excited states. 
In Fig. [3j the increase in the yields of the 2s and 2p 
He + parent ions after the IR pulse as functions of the 
time delay between the two pulses is reported. Both ion 
yields are modulated by substantial quantum beats p3j . 
A similar phenomenon, due to the coherent superposition 
of bound states rather than resonances, was predicted in 
the ionization of C + [13] and has already been observed 
in helium, close to the N=l ionization threshold [25|. In 
the latter case, the authors demonstrated that it is pos- 
sible to control both the timing and the probability of 
ionization. In the present case, on the other hand, we 
show that this holds also for the branching ratio between 
different ionization channels. The ion yields track the 
sp + breathing mode. When fitted with the function 

A • sm(ujt + 0) + Co + c\t + C2^ 2 , 

the two curves in Fig. Ogive the same uo = 0.121(3) a.u., 
which is readily recognized as the energy difference be- 
tween the two lowest sp + resonances, AE = 0.123 a.u. It 
is interesting to note that the two oscillations are out of 
phase by as much as 60°, corresponding to a time delay of 
~ 200 as; hence, the 2s and 2p channels sample different 
stages of the breathing motion of the metastable wave 
packet. Moreover, the population of the spj resonance 
immediately after the pump pulse is larger than that of 



the resonance. Since the spj lifetime is the shortest, 
there is a moment at which the decay rates of the two res- 
onances become comparable. At this point the beating 
between the two resonances is maximal. These features 
suggest the possibility to control the branching ratio of 
the two ionization-excitation channels and, in turn, to 
alter the course of reactions where metastable electronic 
states play a dominant role. 

There are several ways to detect these quantum beats. 
First, the radiative lifetimes of the 2s and the 2p He + 
states differ by several orders of magnitude (~ 1.9 ms [26[ 
and ~ 10 -10 ms (27J respectively), therefore their yields 
can be disentangled by looking at their fluorescence de- 
cay with (2s + 2p) or without (2p alone) the presence of 
an external quenching electric field. Second, since the 
photoelectron angular distribution in the 2s and 2p dif- 
fer, asynchronous beats should be visible in the signal of 
electrons collected along different directions, e.g., along 
the polarization axis and in the plane orthogonal to it. 
The photoelectrons in the N=l and N=2 channels should 
give rise to two well-separated signals already in a veloc- 
ity map imaging spectrometer [28] , because of their very 
different energies. If necessary, the latter could be de- 
tected in coincidence with their excited parent- ion coun- 
terpart by using a reaction microscope |29j . 

In Fig. [2k we show the 2p PDPAD after an IR pulse 
delayed from the probe by 15.53 fs. The probe pulse ion- 
izes the DES generating a short series of peaks above the 
N=2 threshold which interfere with the direct-ionization 
amplitude. In the lapse between the two pulses, the latter 
accumulates a phase which is linear in both the energy 
and the time delay: 

(fdirect(E,t) = <Pdirect(E,t ) + E(t - t )/h. (5) 

As a consequence, prominent interference fringes, with 



characteristic energy spacing AE ~ 27r/(t — to), emerge. 
The sp^ and sp% 1 P° resonances are separated from the 
N = 2 ionization threshold by roughly the energy of three 
and one IR photons, respectively. Since the absorption 
of an odd number of photons by a 1 P° state results in an 
even parity state, the 2p multiphoton ionization ampli- 
tude, created by the IR pulse, should have odd parity just 
above the threshold and change to even parity for photo- 
electron energies around 1.4 eV. Indeed, the interference 
pattern between multiphoton and direct ionization am- 
plitudes (Fig. [2k) is asymmetric with respect to cos close 
to the threshold, and symmetric above E = —0.45 a.u. 
At each increase of the time-delay by half an IR cycle, 
the relative phase between direct and indirect multipho- 
ton amplitudes changes by 7r close to the threshold, while 
it remains the same one photon energy above. Indeed, 
approximately, Fig.2e is the mirror image of Fig.2f. Sim- 
ilar considerations apply for the 2s PDPAD in Figs.[2]>c. 
By measuring the PDPADs at different time delays [25| 
it is possible to recover the ionization amplitude of the 
DES. With additional information on the phase intro- 
duced by the IR field in the ionization amplitude of each 
DES, possibly obtained from simulations, even the orig- 
inal metastable wave packet could in principle be recon- 
structed. 

In conclusion, we have presented evidence that quan- 
tum beating between doubly excited states can be mon- 
itored experimentally, and that it can be exploited with 
the available attosecond pump-probe techniques in order 
to steer the course of atomic photoionization. 
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